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This work presents accurate ab initio determination of the magnetic dipole (Ml) and electric 
quadrupole (E2) hyperfine structure constants for the ground and a few low-lying excited 
states in 67 Zn+ which is one of the interesting systems in fundamental physics. The coupled- 
cluster (CC) theory within the relativistic framework has been used here in this calculations. 
Long standing demands for a relativistic and highly correlated calculations like CC can be 
able to resolve the disagreements among the lifetime estimations reported previously for a few 
low- lying states of Zn + . The role of different electron correlation effects in the determination 
of these quantities are discussed and their contributions are presented. 

I. INTRODUCTION 

The quantum information processing (QIP) is one of the interesting areas in physics which is 
gaining momentum both in theoretical and experimental fronts in the recent years. Mostly, the 
sin gl e vaience ions parting the one, with ^ ground states ate being chosen fot QIP stud.es 
[1[ to encode qubits into the hyperfine levels. These levels are chosen due to their relatively long 
lifetimes against spontaneous decay rates and long phase coherence because of their small energy 
separations. The hyperfine structure studies help us understand the nuclear structure of an atom 
and its influence on the short range wavefunctions correctly Q]. 

The rapid progress in the development of technology involving laser cooling and ion trapping 
has made possible to bring these theoretical ideas to fruition and singly ionized zinc (Zn + ) is one 
of the recent important inclusion in that family [3]. Though there are a few studies of radiative 
lifetimes of Zn + in the literature [4], its hyperfine structures are not studied so far limited to the 
best of our knowledge. Here, we have carried out the magnetic dipole (A) and electric quadrupole 
(B) hyperfine structure studies of 67 Zn + for principle quantum number n=4 states which is wanted 
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for the QIP studies as mentioned above. 

Zn is also one of the important elements in astrophysics, especially for the understanding of the 
post-main sequence evolution of the chemically peculiar stars, in which Zn is either scarce (if not 
non-existing) or over abundant [5J]. The high resolution spectra obtained from GHRS onboard Hub- 
ble Space Telescope has provided vital informations about its abundances [6y. Applications of the 
radiative transitions of this ion in cosmology, stellar dynamics, interstellar medium, nucleosynthesis 



im. 



etc. have been discussed extensively in the literature 

It seems from the reported results that, there exits disagreements in the lifetime estimations 
among the experimental measurements and various theoretical calculations. One distinct feature 
of the lifetime table is the order of the AD fine structure states, i.e., the lifetime of the ADjj 2 state 



should be less than the lifetime of 4L> 3 / 2 state according to both experiments available [9j, |10(]. This 
was not found in any of the ab initio studied so far 9, 11. [12I \\\. Our calculated lifetimes for the 



4D3/2 and 4D 5/ / 2 states which are reported here, have the same order as well as in good agreement 
with the experimental results. 

II. THEORY 

The one-electron reduced matrix elements corresponding to El, Ml and E2 transitions are given 
in these papers [14I . [ijj. The emission transition probabilities (in sec 1 ) for the El, E2 and Ml 
channels from states / to i are given by 



E1 _ 2.0261 x 10 18 El 

A iJ ~ A 3 (2j/ + I) b > l 



m _ 1.11995 x 10 18 E2 



ami _ 2-69735 x 10 13 m 

A ^ ~ A3(2j / + 1) b ' (2 - S) 

where S = f\0\^i)\ 2 is the transition strength for the operator O (in a.u.) and A (in(A)) is 
the corresponding transition wavelength. The lifetime of a particular state is the reciprocal of the 
total transition probability arising from all possible spontaneous electromagnetic transitions from 
the state to all the lower energy levels. 



x- (2 - 4) 
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The interaction between the electromagnetic multipole moments of the electrons and the elec- 
tromagnetic field created at the site of the nucleus is termed as hyperfine interaction and the 
corresponding Hamiltonian is given by 1|| 



i^^MW-TM (2.5) 

k 

where and are the spherical tensor operators of rank k in the nuclear and electronic 

spaces, respectively. The k=l and 2 terms of the expansion represent the magnetic dipole and 
electric quadrupole interactions, respectively. 

The diagonal hyperfine interaction constants can be written as 

/ 



a = %{,jjW\uj) = % 



J 1 J 

■JO J 



(JHT^HJ), (2.6) 



and 



B = 2Q( 7 JJ\T^\ 7 JJ) = 2Q \ J 2 J \ (J||T( 2 )||J), (2.7) 

I -J J J 

where /, J are the total angular momentums of nucleus and electrons; /i/ and Q are magnetic 
dipole and electric quadrupole moments of the nucleus, respectively. The and T^ 2 ) operators 
are defined as 

T (i) = £ -ieV2rr 2 ai ■ Cfo\) (2.8) 

i 

and 

T (2) = J2-err*C 2 q (n), (2.9) 

i 

where, C q = ^/47r/(2/c + l)Yk q with Y^ q being the spherical harmonic functions. 

In the first-order perturbation theory, the hyperfine energy Ehf s {J) of the fine-structure state 
\JMj) is the expectation value of the corresponding hyperfine interaction Hamiltonians in that 
state. The energies corresponding to the magnetic dipole and electric quadrupole hyperfine tran- 
sition are defined as 



E M i = AFT/2, 



(2.10) 
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and 



E. 



Q2 



B 3K(K + 1) - 4/(/ + I) J (J + 1) 



(2.11) 



2 2/(2/ - 1)2J(2J - 1) 
where K = 2(1 ■ J) = F(F + 1) - 1(1 + 1) - J( J + 1) with F = I + J. Here we have neglected 
higher order hyperfine interactions. 

The basic formalism of the valence universal coupled-cluster (CC) method was developed more 



than two decades before 
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however a suitably relativistic version of this approach has 
been successfully employed to obtain the various properties accurately in different single valence 
atomic systems only recently 21 



251 ] . Here we just outline the method applied in this 



calculation of the wavefunctions of Zn + accurately. 

The single valence CC theory extended for the relativistic framework and is based on the no- 

n 

virtual-pair approximation with Dirac-Fock orbitals .18]. The concept of the common vacuum 
for both the closed-shell N and open-shell N + 1 electron systems allows to formulate a direct 
method of excitation energies. The dynamical electron correlation effects are introduced through 
the valence-universal wave-operator £l v 
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la ] for the state with v as the valence orbital is written 



in the normal ordered form as, 



e r {e s *}, 



(2.12) 



where cluster operator T represents excitations from the occupied core orbitals of the closed shell 
system Zn ++ and S represent the core-valence and valence-valence excitations. Dominant among 
these correlations are pair correlations and core polarizations. The Dirac-Coulomb Hamiltonian 
dressed with the excitation cluster operators T and S v are then diagonalized within the model space 
constructed from the core and valence orbitals to obtain the desired eigenvalues and eigenvectors 
[20I ] . In this work, a leading order triple excitations are included in the open shell CC amplitude 
evaluation by an approximation that is similar in spirit to CCSD(T) method 261 ] . 
The expectation value of any operator O can be expressed in the CC method as 

<*„|0|*„> 



o 



($ v \{l + Sj}e T *Oe T {l + S v }\$ v ) 
(0> v \{l + Sj}e Ti e T {l + S v }\* v ) 



(2.13) 



III. RESULTS AND DISCUSSIONS 



We have used Gaussian- type orbitals (GTO) to calculate the DF wavefunctions \$df) a s given 



111 



271 ] using the basis functions of the form [2J 
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Gi,*(r) = r k >e- a * r (3.1) 

where k = 0, 1, 2, 3, .... for s, p, d, f type orbital symmetries respectively. The large and small 

components of the relativistic GTOs satisfy the kinetic balance condition [31]. The exponents are 
determined by the even tempering condition; i.e., for each symmetry exponents are assigned as 

an = a p i ~ 1 * = 1,2, N (3.2) 

where N is the number of basis functions for the specific symmetry. In this calculation, we have 
used «o = 0.00831 and (3 = 2.99. The number of basis functions used in the present calculation 
are 32, 32, 30, 25, 20 for I = 0, 1, 2, 3, 4 symmetries, respectively. 

TABLE I: Radiative lifetimes (ns) for different low-lying states of 67 Zn + . 

State Experiment Other theories This work 

4 2 P 1/2 2.6(3) [9], 3.05(4) [35], 2.1(4) [36], (2.19, 2.01) [9], 2.151 [11], (1.97, 2.38) [13], 2.43 

2.5(2) [4], 2.7 [37], 2.524 [12], 2.2 [36], 3.2 [38], 2.41 [39] 

4 2 P 3/2 2.6(3) [9], 2.07(0.2) [40], 2.1(3) [36], (2.1, 1.91) [9], 2.036 [U], (1.85, 2.23) [13], 2.30 

3.1(4) [35], 3.0(3) [10], 2.5(2) [4], 2.5 [37], 2.386 [12], 2.0 [36], 2.9 [38], 2.27 [39] 

S 2 ^^ 1.7(14) [9], 2.0(2) [41], 1.8(2) [40], (1.86, 1.99) [9], 2.5 [36], 2.468 [11], 2.08 

3.85(7) [35] 2.07 [13], 1.4 [12] 

4 2 L> 3/2 1.8(2) [9J, 1.8(4) [10] (1.27, 1.16) [9], 2.262 [11], 1.363 [12], 1.39 [13] 1.31 

4 2 D 5/2 1.44(12) [9J, 1.40(15) [10] (1.31, 1.21) [9], 1.285 [11], 1.388 [12], 1.33 [13] 1.04 

We report our calculated lifetime results along with the available calculated and measured 
results in Table I. It is apparent from the table that, there are large disagreements among the 
earlier results and also the measurements are not very precise. In our previous work (here after 
referred to as paper I [32]), we have presented ionization energies, allowed and forbidden transition 
amplitudes of the same system considered in the present work and their astrophysical applications 
are emphasized Our results in paper I are in excellent agreement with the experimental 

measurements. With this spirit we have computed the lifetime calculations in the present work 
which show moderately good agreement with the available experimental results. There are two 
recent experiments for the lifetime estimations of the fine structure states of 4Z) [l(J and both 
show that, the lifetime of the 4L> 5 / 2 state is shorter than 4Z? 3( / 2 state. However, many of the earlier 
theoretical calculations in the literature show the opposite trend, whereas, our CCSD(T) results 
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FIG. 1: Decay channels for the first few low- lying excited states of Zn + . The lines of different types 
correspond to different electromagnetic (multipole) transitions. 

not only show the same trend as that of the experimental results, but also the ratio of their lifetimes 



of lifetime, as seen in Fig. 1. 

The computed values of the magnetic dipole hyperfine structure constants (A) for the ground 
state and a few low-lying excited states of 67 Zn + are given in Table II. Neither the calculations 
nor the measurements of the hyperfine constants A for all the states, except for the ground 
state, considered here are available in the literature known to our knowledge. The important 
many-body correlation contributions to the total A^ values are included in our work through 
relativistic CC theory. We have used fij = 0.87547 and / = 5/2 [33] in our calculation. From the 
differences of the DF and the total CC results it is evident that the electron correlation effects 
in the calculation of A are quite large and vary from (2-430)% among the different low- lying 
states. The core correlation effects are significant inthe 4p 1( / 2 , 4d 3 / 2 and 4d 5 / 2 states; especially 
for the last two states these effects are larger than the DF contributions. The lowest order pair 
correlation and core polarization effects tabulated here highlight their important contributions 
which are comparable to the DF contributions, especially to note is the cancellation effect of core 
polarization in the case of 4p 3/2 state. The large effects of S\ v OS 2v (4.93 MHz for 4d 5/2 and 40.0 
MHz for 4s state) are observed in these cases. We have used the expression (2.10) to calculate 
the ground state hyperfine energy separation which turns out to be 7018.743 MHz. Panigrahy 
et al. also have calculated the same using relativistic linked-cluster many body perturbation 
theory and get 7.2 GHz [3J] which is in good agreement with our result. The hyperfine energy 
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TABLE II: Magnetic dipole hyperfine constant (A) of different low-lying states of 67 Zn + in MHz. 



otate 


1 J 1 


Core Correlation 


Pair Correlation 


Core Polarization 


Norm 


loiai 


4Si/ 2 


1835.57 


32.64 


345.11 


100.68 


-53.27 


2339.58 


4A/2 


271.77 


146.90 


65.33 


37.70 


-9.65 


526.68 


4^3/2 


49.94 


2.07 


11.68 


-15.84 


-0.92 


50.74 


4£>3/2 


6.33 


6.46 


0.22 


1.71 


-0.07 


15.57 


4^5/2 


2.72 


4.54 


0.72 


1.57 


-0.07 


14.56 



TABLE III: Electric quadrupole hyperfine constant (B) of different low-lying states of 67 Zn + in MHz. 



State 


DF 


Core Correlation 


Pair Correlation 


Core Polarization 


Norm 


Total 


4F 3 /2 


40.01 


0.72 


9.43 


10.74 


-1.12 


61.74 


4£>3/2 


1.68 


0.23 


0.37 


2.71 


-0.02 


5.04 


4^5/2 


2.40 


-0.22 


0.48 


3.15 


-0.03 


7.35 



separation lies in the microwave region of the electromagnetic spectrum, which suggest that 67 Zn + 
can be proposed as the new frequency standard in microwave region, however it needs further 
investigation about its stability and accuracy of estimation etc. 

The computed values of the electric quadrupole hyperfine structure constants (B) for a few 
low-lying excited states are given in_Table III. In this calculation the electric quadrupole moment 
of the nucleus, Q = 0.150 is used 



331 ]. It may be noted that the effects of pair correlation and 
core polarization effects are stronger than the core correlation effects; in particular, the core 
polarization effects for the D states are stronger than the DF contributions. 



IV. CONCLUSION 



In this work, we have determined the hyperfine structure constants A and B of the ground state 
and a few low- lying excited states in 67 Zn + using the relativistic coupled-cluster theory. We have 
also calculated the hyperfine energy separation for the ground state, which is 7018.743 MHz. There 
is no experimental result available for the hyperfine energy separation for the ground state, which 
seems to be an important candidate for QIP studies. Also, 67 Zn + can be considered as one of the 
promising candidates for the frequency standard in the microwave region. We have also determined 
the lifetimes of the low-lying state sin Zn + , which are in good agreement with experimental results. 
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Especially our calcu 
in the experiments 



ated lifetimes of AD fine structure states explain the same trend as observed 



<], fl i.e., 



the lifetime of 4D 5/ / 2 state is shorter than the 4Z) 3 / 2 state, unlike 
many other theoretical results which show opposite trend and also the ratio of their lifetimes in our 
calculation is in excellent agreement with the experimental result. This suggests that the relativistic 
CC method applied in the present work and our numerical approach in obtaining the wavefunctions 
of the system considered are more accurate and reliable. 
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